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ABSTRACT 

We present high-redshift predictions of the star formation rate distribution func¬ 
tion (SFRDF), UV luminosity function (UVLF), galactic stellar mass function 
(GSMF), and specific star-formation rates (sSFRs) of galaxies from the latest ver¬ 
sion of the Munich semi-analytic model L-Galaxies. We find a good fit to both 
the shape and normalisation of the SFRDF at z = 4 — 7, apart from a slight under¬ 
prediction at the low SFR end at z = 4. Likewise, we find a good fit to the faint number 
counts for the observed UV LF; at brighter magnitudes our predictions lie below the 
observations, increasingly so at higher redshifts. At all redshifts and magnitudes, the 
raw (unattenuated) number counts for the UVLF lie above the observations. Because 
of the good agreement with the SFR we interpret our under-prediction as an over¬ 
estimate of the amount of dust in the model for the brightest galaxies, especially at 
high-redshift. While the shape of our GSMF matches that of the observations, we lie 
between (conflicting) observations at z = 4 — 5, and under-predict at z = 6 — 7. The 
sSFRs of our model galaxies show the observed trend of increasing normalisation with 
redshift, but do not reproduce the observed mass dependence.Overall, we conclude 
that the latest version of L-Galaxies, which is tuned to match observations at z ^ 3, 
does a fair job of reproducing the observed properties of galaxies at z ^ 4. More work 
needs to be done on understanding observational bias at high-redshift, and upon the 
dust model, before strong conclusions can be drawn on how to interpret remaining 
discrepancies between the model and observations. 

Key words: galaxies: evolution - galaxies: formation - galaxies: luminosity function, 
mass function - galaxies: high-redshift - ultraviolet: galaxies 


1 INTRODUCTION 

With the installation of Wide Field Camera 3 (WFC3) 
on the Hubble Space Telescope (HST) in 2009 it is now 
possible to identify statistically useful and robust sam¬ 
ples of star forming galaxies in the early Universe (z > 4, 
Oesch et al. 2010; Bouwens et al. 2010a; Bunker et al. 2010; 
Wilkins et al. 2010; Finkelstein et al. 2010; Mc.Lure et al. 
2010; Wilkins et al. 2011a; Lorenzoni et al. 2011; 
Bouwens et al. 2011; McLure et al. 2011; Finkelstein et al. 
2012a; Lorenzoni et al. 2013; McLure et al. 2013; 
Duncan et al. 2014; Finkelstein et al. 2014). In recent 
years a tremendous effort has been dedicated to quantifying 
the photometric and physical properties, such as star 
formation rates and stellar masses, of these galaxies. As we 
continue to dig deeper, with the first sources now identified 
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at z « 10 (e.g. Oesch et al. 2012; Ellis et al. 2013), and with 
the launch of the James Webb Space Telescope (JWST) in 
the next few years, we will further be able to constrain the 
physics of galaxy formation and evolution in this critical 
epoch of the Universe’s history. 

Although it lasts less than 0.8 Gyr, the period of the 
Universe between z = 7 and z = 4 is important to study 
because it defines an epoch of interesting galaxy formation 
and evolution activity. The start of this period marks the 
end of the epoch of reionization; galaxies are starting to 
ramp up their metal and dust production; and we are find¬ 
ing evidence of the first quasars. While astronomy is unique 
in allowing us to observe the Universe at these early times, 
theoretical modelling is required to interpret those observa¬ 
tions in terms of an evolving galaxy population. The rapidly 
advancing observational constraints on the physical proper¬ 
ties of galaxies in the early Universe provides an opportunity 
to further test and refine these galaxy formation models. 
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The most well studied property of the galaxy population 
at high-redshift (in part due to its accessibility) is the rest- 
frame ultraviolet (UV) luminosity function (LF). Because 
of the link between the UV luminosity of galaxies and their 
star-formation rates, the observed UV LF provides an impor¬ 
tant constraint on star-formation activity in the early Uni¬ 
verse. While early observational results were based on only 
small samples (Bouwens et al. 2008; Bouwens et al. 2010b, a; 
Bunker et al. 2010; Oesch et al. 2009; Oesch et al. 2010; 
Ouchi et al. 2009; Wilkins et al. 2011b; Robertson et al. 
2010; Dunlop et al. 2010; Lorenzoni et al. 2011), we have 
slowly begun building larger catalogues, first with 200 — 
500 galaxies (Finkelstein et al. 2010; Bouwens et al. 2011; 
McLure et al. 2013), with the most recent observations hav¬ 
ing almost 1000 galaxies at z ^ 7 (Bouwens et al. 2015; 
Finkelstein et al. 2014). 

While the intrinsic UV luminosity is known to be a use¬ 
ful diagnostic of star-formation activity (e.g. Wilkins et al. 
2012), it is susceptible to even small amounts of dust (Auv ~ 
10 x E(B — V)). Direct comparison of the observed UV lu¬ 
minosity function with predictions from galaxy formation 
models is then sensitive to the reliability of the dust model 
(which has to account for the creation and destruction of 
dust, and its effect on the intrinsic spectral energy distribu¬ 
tion). 

Whilst challenging, it is observationally possible to con¬ 
strain the dust obscuration and thus determine the true (or 
intrinsic) star formation activity, even in distant galaxies. 
Starlight that is absorbed by dust is reprocessed and emitted 
in the rest-frame mid/far-IR. Combining the star-formation 
rate inferred from the observed UV with that inferred from 
the mid/far-IR. omission then provides a robust constraint 
on the total (or intrinsic) star-formation activity. Observa¬ 
tional constraints on the rest-frame mid/far-IR emission in 
high-redshift galaxies are, however, challenging due to the 
significantly lower flux sensitivity and poorer spatial reso¬ 
lution of facilities operating at these wavelengths. Thus far 
there is only a single galaxy individually detected in the far- 
IR at z > 6 (Riechers et al. 2013). This is, however, likely to 
rapidly improve with the completion of the Atacama Large 
Millimetre Array (ALMA). 

One alternative to using far-IR/sub-mm observations 
is to take advantage of the relationship between the 
rest-frame UV continuum slope /3, which is easily accessible 
even at z ~ 10 (Wilkins et al. submitted) and the UV 
attenuation (first applied by Meurer, Heckman & Calzetti 
1999). The measurement of (5 in high-redshift galaxies 
has, in recent years, been the focus of intense study 
(e.g. Stanway, McMahon & Bunker 2005; Bouwens et al. 
2009; Bunker et al. 2010; Bouwens et al. 2010b; 
Wilkins et al. 2011b; Dunlop et al. 2012; Bouwens et al. 
2012; Finkelstein et al. 2012b; Castellano et al. 2012; 
Rogers, McLure & Dunlop 2013; Wilkins et al. 2013; 
Bouwens et al. 2014). Measurements of the UV continuum 
slope have been used to effectively correct the observed UV 
luminosity function and thus determine the star-formation- 
rate distribution function (e.g. Smit et al. 2012). It is 
important to note, however, that this relation is sensitive 
to a number of assumptions (see Wilkins et al. 2012, 2013) 
which introduce both systematic biases and increase the 
scatter in individual observations. 

By combining space (from Hubble) and ground-based 


near-IR observations (< 2/xm) from the Infra-red Array 
Camera (IRAC) aboard the Spitzer Space Telescope it is 
possible to probe the rest-frame to optical spectral energy 
distributions (SEDs) of galaxies at high redsliift. This is crit¬ 
ical to deriving robust stellar masses and thus the galaxy 
stellar mass function (GSMF). The measurement of stellar 
masses at high-redshift is, unfortunately, affected by vari¬ 
ous issues, including: the low sensitivity of the IRAC ob¬ 
servations; assumptions regarding the star formation and 
metal enrichment history of these galaxies; and the effects 
of strong nebular emission (e.g. Wilkins et al. 2013). De¬ 
spite these obstacles, several groups have now attempted to 
measure the galaxy SMF in the high-redshift Universe (e.g. 
Stark et al. 2009; Labbe et al. 2010; Gonzalez et al. 2011; 
Yan et al. 2012; Duncan et al. 2014) permitting a direct 
comparison with galaxy formation models. 

The Munich semi-analytic model of galaxy formation 
(latest version Henriques et al. 2014), also known as L- 
Galaxies, has had a lot of success over the past decade in 
predicting various properties of galaxies, such as the stellar- 
mass and luminosity functions both in the local Universe 
and out to redshift z = 3 (Henriques et al. 2013). In this 
paper we extend these predictions to higher redshift with¬ 
out altering any of the model parameters (except to modify 
the redshift-dependence of the dust model, as described in 
§2.2.1 below). In that sense, the results presented here may 
be considered predictions of the model. 

This paper is structured as follows: in Section 2 we de¬ 
scribe the relevant parts of our semi-analytical model, high¬ 
lighting the changes in the latest version; in Section 3 we dis¬ 
cuss our high redshift predictions for the star-formation-rate 
distribution function and UV luminosity function, followed 
by the stellar mass function in Section 4. In Section 5 we dis¬ 
cuss the relationship between the specific star formation rate 
and the stellar mass and in Section 6 we give a brief overview 
of the evolution of these properties at high-redshift. We 
conclude our work in Section 7. Throughout this paper we 
adopt a Chabrier initial mass function (Chabrier 2003) and 
use the latest Planck cosmology (Planck Collaboration et al. 
2014). Number densities are presented per co-moving vol¬ 
ume, (h _1 Mpc) 3 . 

This paper makes use of detailed predictions from 
the new model of L- Galaxies outlined in Henriques et al. 
(2014), which have been made publicly available. 1 The 
binned data used to make the plots in this paper have also 
been made available online. 2 

2 THE MODEL 

2.1 L- Galaxies 

Semi-analytic models (SAMs) provide a relatively inexpen¬ 
sive method of self-consistently evolving the baryonic com¬ 
ponents associated with dark matter merger trees, derived 
from V-body simulations or Press-Schechter calculations. 
The term semi-analytic comes from the use of coupled dif¬ 
ferential equations (rather than numerical calculations), to 


1 http://gavo.mpa-garching.mpg.de/MyMillcnnium/ 

2 http://astronomy.sussex.ac.uk/~sc558/ 
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Figure 1. Predicted stellar mass functions at redshift 2 f5 0 (top left); z ~ 1 (top right); z zz 2 (lower left) and 21 3 (lower right). Solid 

black lines show the stellar mass functions predicted by our model. This figure is reproduced from Henriques et al. (2014) Figures 2 and 
Al, and we direct the reader there for a more detailed description. We include it here to highlight how well the model works at lower 
redshifts in predicting key observables such as the SMF. Observations are taken from several surveys; SDSS (Baldry, Glazebrook & Driver 
2008; Li & White 2009) and GAMA (Baldry et al. 2012) at z=0; and Marchesini et al. (2009), Spitzer-COSMOS (Ilbert et al. 2010), 
NEWFIRM (Marchesini et al. 2010), COSMOS (Dominguez Sanchez et al. 2011), ULTRAVISTA (Muzzin et al. 2013; Ilbert et al. 2013) 
and ZFOURGE (Tomczak et al. 2014) at higher redshifts. 


follow the evolution of galaxy formation physics determin¬ 
ing the properties of gas and stars. Physics commonly found 
in most SAMs include descriptions of: (1) primordial infall 
and the impact of an ionizing UV background; (2) Radia¬ 
tive cooling of the gas; (3) Star Formation recipes; (4) Metal 
enrichment; (5) Super-massive black hole growth; (6) Super¬ 
novae and AGN feedback processes; (7) The impact of en¬ 
vironment and mergers including galaxy morphologies and 
quenching. 

The Munich SAM, or L-Galaxies, (Springel et al. 
2001; De Lucia, Kauffmann & White 2004; Springel et al. 
2005; Croton et al. 2006; De Lucia & Blaizot 2007; 
Guo et al. 2011; Guo & White 2013; Henriques et al. 
2013) has been developed over the years to include most 
of the relevant processes that affect galaxy evolution. In 


this work we use its latest version, Henriques et al. (2014), 
and direct the reader to the appendix of that paper for a 
detailed description of the model. Of most relevance to this 
paper are the adoption of the Planck year 1 cosmology 
and a modified gas-to-dust relation, partly motivated by 
the work presented in this paper. The model parameters 
were constrained using the abundance and passive fractions 
of galaxies at z ^ 3; and the model has successfully 
reproduced key observables at these redshifts, such as the 
luminosity and stellar mass functions. We highlight this 
fact in Figure 1, which is a reproduction of the SMF at 
z € {0,1,2,3} from Henriques et al. (2014). We direct the 
reader to Figures 2 and Al and the related text of that 
paper for a more detailed discussion, but we highlight how 
well the model can explain the observed evolution in the 
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SMF at these redshifts, over the mass range constrained by 
observers. 


2.2 Dust Extinction Model 

Actively star-forming galaxies are known to be rich in dust. 
This can have a dramatic effect on their emitted spectrum 
since dust significantly absorbs optical/UV light while hav¬ 
ing a much milder effect at longer wavelengths. As a re¬ 
sult, dust-dominated galaxies will generally have red colours 
even if they are strongly star-forming. For that reason, we 
summarise the dust model of Henriques et al. (2014) here: 
a fuller description can be found in Section 1.14 of the sup¬ 
plementary material in that paper. 

We considering dust extinction separately for the diffuse 
interstellar medium (ISM) and for the molecular birth clouds 
(BC) within which stars form. The optical depth of dust as 
a function of wavelength is computed separately for each 
component and then combined as described below. We do 
not at present attempt to compute the detailed properties 
of the dust particles or the re-emission of the absorbed light. 


2.2.1 Extinction by the ISM 


The optical depth of diffuse dust in galactic disks is assumed 
to vary with wavelength as 


where 


ISM 

T\ 


= (l + z)- 1 




(Nh) 


2.1 x 10 21 atoms cm -2 


(Nh) 


A/cold 

1.4 m.p7r(nT?g aSi d )^ 


( 1 ) 

( 2 ) 


is the mean column density of hydrogen. Here -R gas ,d is the 
cold gas disk scale-length, 1.4 accounts for the presence of 
helium and a = 1.68 in order for (Nh) to represent the mass- 
weighted average column density of an exponential disk. Fol¬ 
lowing the results in Guiderdoni (1987), the extinction curve 
in eq. (1) depends on the gas metallicity and is based on 
an interpolation between the Solar Neighbourhood and the 
Large and Small Magellanic Clouds: s = 1.35 for A < 2000 
A and s = 1.6 for A > 2000 A. The extinction curve for 
solar metallicity, (A\/Ay)z @ , is taken from Mathis (1983). 

The redshift dependence in eq. (1) is significantly 
stronger than in previous versions of our model ((1+z) -0 ' 5 in 
Kitzbichler & White (2007) and (1 + z)~ 0 ' 4 in Guo & White 
(2009)). The dependence implies that for the same amount 
of cold gas and the same metal abundance, there is less dust 
at high redshift. The motivation comes both from observa¬ 
tions (Steidel et al. 2004; Quadri et al. 2008) and from the 
assumption that dust is produced by relatively long-lived 
stars. However, it may also be that this redshift dependence 
has to be introduced as a phenomenological compensation 
for the excessively early build-up of the metal content in 
model galaxies. In practice it has been included simply to 
give an approximate match to the low extinctions of high- 
redshift galaxies as inferred from their observed UV slopes 
(Bouwens et al. 2012), and to the UV luminosity function, 
as described below. 


2.2.2 Extinction by molecular birth clouds 

This second source of extinction affects only young stars that 
are still embedded in their molecular birth clouds, for which 
we assume a lifetime of 10 Myr. The relevant optical depth 
is taken to be 

where /r is given by a random Gaussian deviate with mean 
0.3 and standard deviation 0.2, truncated at 0.1 and 1. 


2.2.3 Overall extinction curve 


In order to get the final overall extinction, every galaxy is 
assigned an inclination, 8, given by the angle between the 
disk angular momentum and the 2 -direction of the simula¬ 
tion box, and a “slab” geometry is assumed for the dust in 
the diffuse ISM. For sources that are uniformly distributed 
within the disk then the mean absorption coefficient is 


Ai SM = —2.5 lo gl 


„ __ISM c 

1 — exp A 


r ISM 


sec# 


(4) 


Emission from young stars embedded within birth clouds is 
subject to an additional extinction of 

= -2.51og 10 (exp _r ^ C ) . (5) 

The standard L- Galaxies output does not attempt 
to model the attenuation of light by the intergalactic 
medium. However, this is done in post-processing for 
the lightcones published in the Millennium Run Observa¬ 
tory 3 (Overzier et al. 2013). In this paper, however, we ne¬ 
glect intergalactic attenuation. 


3 RECENT STAR FORMATION 

In this section we investigate the star formation rate, and the 
related UV luminosity function, at redshifts z £ {4, 5,6, 7}. 

Figure 2 shows the star-formation-rate distribution 
function (SFRDF) as predicted by our model alongside mea¬ 
surements from Smit et al. (2012) (converted to our fiducial 
Chabrier IMF) and Duncan et al. (2014). 

Comparing with the Smit et al. (2012) measurements 
at redshifts z sj 5 — 7 we find generally good agreement. At 
these redshifts, the Duncan et al. (2014) measurements are 
generally higher than both the model and the Smit et al. 
(2012) results. This is particularly true for the most mas¬ 
sive galaxies, though we note that the quoted observational 
uncertainties can be very large. 

At z ~ 4, however, our model under-predicts the num¬ 
ber of galaxies for log 10 (SFR//C 2 MQ) < 1 when compared 
to both sets of observations (which are consistent with one 
another at this redshift). The cause of the discrepancy is 
unclear, though may be a consequence of our model under¬ 
estimating the contribution to the SFR from merger-driven 
activity. 

3 The Millennium Run Observatory, or MRObs, allows you to 
observe our semi-analytic galaxy formation model through the 
use of ‘virtual telescopes’. 
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Figure 2. Predicted star formation rate distribution functions at redshift z ss 4 (top left); z ~ 5 (top right); 2 fa 6 (lower left) and 2 fa 7 
(lower right). In each instance we use the closest available snapshot from our L-Galaxies run of z=3.95, 5.03, 5.82 and 6.97 respectively. 
Solid black lines show the star formation rate distribution function predicted by our model. Our 2 = 4 star formation rate distribution 
function is repeated at higher redshifts as a grey dot-dash line for comparison. Observations are taken from Smit et al. (2012), converted 
to a Chabrier IMF, and Duncan et al. (2014). 


3.1 The UV Luminosity Function 

We present the UV luminosity function predicted by our 
model in Figure 3 alongside recent observational estima¬ 
tions at high-redshift (Bouwens et al. 2015; Duncan et al. 
2014; Finkelstein et al. 2014; Bowler et al. 2014a, b). The 
solid black line shows our prediction for the attenuated UV 
luminosity function; the attenuated UV LF at 2 = 4 is also 
shown on subsequent plots for comparison. The dashed line 
shows our intrinsic UV luminosity function, with no dust 
model being applied. 

We find a good fit to the faint number counts: Muv > 
—20 for 2 = 4 — 6 and Muv > —19 for 2 = 7. At brighter ab¬ 
solute magnitudes, the model counts fall below the observed 
ones. Note, however, that the raw counts, before dust at¬ 
tenuation, lie above the observations. Given that we saw a 
good fit in Section 3 between predicted and observed SFRs, 


then this points to a difference in the dust model between 
the two. 

To better understand this, we quantify in Figure 4 the 
attenuation required (as a function of the intrinsic UV abso¬ 
lute magnitude) to reconcile the raw L-Galaxies data with 
observations. We do this by comparing observed, M$ j0 b a , 
and intrinsic, M$,i nt , absolute magnitudes below which we 
achieve a particular cumulative number density, 4? of galax¬ 
ies: 



where 4> is the usual differential number density of galaxies. 
The attenuation is then Auv = M<j j0 b s — M$,i n t. 

The dust attenuation required to match the observa¬ 
tion (as a function of the intrinsic absolute magnitude) is 
shown in Figure 4. The black, solid line shows the attenua- 
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Figure 3. Predicted rest—frame (1500 A) UV luminosity functions at redshift z ~ 4 (top left); z ~ 5 (top right); z « 6 (lower left) and 
z ~ 7 (lower right). In each instance we use the closest snapshot available to use from our L-Galaxies run of z=3.95, 5.03, 5.82 and 
6.97 respectively. Solid black lines shows the L-Galaxies prediction for the attenuated UV luminosity function using the dust extinction 
model outlined in §2.2.1. The dashed black line is the L-Galaxies prediction of the intrinsic UV luminosity function, with no dust model 
applied. Our z = 4 attenuated UV luminosity function is repeated at higher redshifts as a grey dot-dash line for comparison. Observations 
are taken from Bouwens et al. (2015), Duncan et al. (2014) and Finkelstein et al. (2014), and at high mass from Bowler et al. (2014a) 
(z = 6) and Bowler et al. (2014b) (z = 7). 


tion currently implemented in L-Galaxies, as described in 
Section 2.2. 

As expected, the built-in attenuation matches that from 
the Duncan et al. (2014) data fairly well. The other data 
sets show a shallower slope: the attenuation is reasonable, 
perhaps even under-estimated in the faintest galaxies, but 
is strongly over-estimated in the brightest galaxies and in¬ 
creasingly so at high redshift. 

It is important to stress that while we are presenting 
the results for all the objects within our simulation, obser¬ 
vational samples (such as those employed by Bouwens et al. 
2015; Duncan et al. 2014; Finkelstein et al. 2014), are bi¬ 
ased and may not truly capture the full galaxy population 
at these redshifts. Indeed, a defining characteristic of the 
Lyman break technique, which is regularly used to iden¬ 


tify galaxies in the high redshift universe, is that it prefer¬ 
entially selects blue rest-frame UV bright sources, i.e star 
forming galaxies with low UV dust attenuation (Auv < 2). 
Very dusty galaxies, or those with little to no star forma¬ 
tion would then be missed in typical Lyman break galaxy 
searches (e.g. HFLS3, a very dusty intensely star forming 
galaxy at z « 6.3 Riechers et al. 2013). The degree to which 
this is a concern at high-redshift is difficult to assess, largely 
due to the lack of sensitive far-IR and sub-mm imaging 
which is critical to identify heavily obscured systems. 

Given the current observational uncertainties, we con¬ 
clude that the simple, empirical dust extinction model cur¬ 
rently built into L-Galaxies does a reasonable job, although 
it could be refined to match particular data sets if required. 
In the future, we intend to implement a more physically- 
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Figure 4. This figure shows the amount of dust attenuation required to move our intrinsic UV luminosity function (the dashed, black lines 
in Figure 3) to match different observational data sets (Bouwens et al. 2015; Duncan et al. 2014; Finkelstein et al. 2014), as a function 
of unattenuated absolute UV magnitude. The solid, black line shows the attenuation built into the L-Galaxies model as described in 
Section 2.2. 


motivated dust model: we note that the current model has 
prompt recycling, and this could be an issue at these early 
times when the age of the Universe is just 1.5 Gyr at 2 = 4 
and less than 1 Gyr for z > 6. A delayed chemical en¬ 
richment model has been implemented in L-Galaxies by 
Yates et al. (2013) and we intend to incorporate that into 
the Henriques et al. (2014) model in future work. 


4 GALAXY STELLAR MASS FUNCTION 

We present the Galaxy Stellar Mass Function (GSMF) at 
2 € {4, 5,6, 7} predicted by our model in Figure 5 along¬ 
side recent observational estimates at high-redshift from 
Gonzalez et al. (2011) and Duncan et al. (2014). 

It is important to first note that the observationally- 
derived mass functions presented in Figure 5 are inconsis¬ 
tent with each other at z ~ 4 — 5. One possible source 
(see Duncan et al. (2014) for a wider discussion) of this 


discrepancy is the effect of nebular emission which was in¬ 
cluded in Duncan et al. (2014) but not in Gonzalez et al. 
(2011). Galaxies in the high-redshift Universe are expected 
(Wilkins et al. 2013) and inferred (e.g. Smit et al. 2014) to 
exhibit strong nebular emission which can strongly affect 
the measured stellar mass-to-light ratios and thus masses 
(Wilkins et al. 2013). The accuracy/precision of stellar mass 
estimates are also affected by the lower sensitivity and an¬ 
gular resolution of the Spitzer /IRAC imaging. 

Given the above observational uncertainties, it is grat¬ 
ifying that the model predictions split the two observa¬ 
tional measurements at 2 = 4. There is a hint that the 
change in slope at the “knee” of the mass-function (M^ nee ~ 
3 x 10 9 h -2 M@) may be sharper in the model than the obser¬ 
vations, but the observational error bars are growing at this 
point and so it is hard to draw firm conclusions. As we move 
to higher redshifts, however, the model predictions and the 
observations gradually diverge as follows: (i) the normalisa- 
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Figure 5. Predicted stellar mass functions at redshift z ~ 4 (top left); z ~ 5 (top right); z ~ 6 (lower left) and z fs 7 (lower right). In 
each instance we use the closest snapshot available to use from our L-Galaxies run of z=3.95, 5.03, 5.82 and 6.97 respectively. Solid 
black lines show the stellar mass functions predicted by our model. To indicate the possible effect of uncertainties in the observational 
stellar mass determinations, we also show as a red dot-dasli line the stellar mass function convolved with a gaussian of standard deviation 
0.3 dex. Our z = 4 stellar mass function is repeated at higher redshifts as a grey dot-dash line for comparison. Observations are taken 
from Gonzalez et al. (2011), converted to a Chabrier IMF, and Duncan et al. (2014). 


tion at Mknee declines more rapidly with increasing redshift 
in the models than in the observations; (ii) the slope of the 
mass function above the knee is steeper in the models than 
in the observations. 

The exact cause of these discrepancies is difficult to as¬ 
sess. One possibility is that it reflects a deficiency in the 
model; on the other hand it may reflect a systematic bias 
in the observations. This has been discussed at low redshift 
(z = 0 — 3) in Appendix C of Henriques et al. (2013). It 
seems probable that the uncertainties on the individual stel¬ 
lar masses could have been underestimated, and that can 
strongly boost the inferred number of galaxies in regions 
where where the mass function is particularly steep. As an 
example of the possible magnitude of this effect, we show in 
Figure 5 the result of convolving with a gaussian of standard 
deviation 0.3 dex, similar to that required at low redshift. 


This largely reconciles the observed and predicted slopes of 
the mass function, but the normalisation remains too low 
at z = 7. Understanding the source of this discrepancy is 
a focus of an additional work in progress (Wilkins et al. 
in-prep). 

Recent hydrodynamic simulations, particularly Illus- 
tris (Vogelsberger et al. 2014) and Eagle (Schaye et al. 
2015), have begun making predictions of observables at 
high redshift (Genel et al. 2014; Furlong et al. 2014). Like 
L-Galaxies, both Illustris and Eagle make predictions 
at high redshift by only using observational constraints at 
lower redshift. Both simulations are similar to ours in the 
prediction of the GSMF at z = 6 — 7 in that we all un¬ 
der predict the abundance of high mass galaxies (> 1O 9 M 0 ) 
at these redshifts, although Eagle better match the obser¬ 
vations at a = 5 across the entire mass range. Whilst both 
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Figure 6. Predicted specific star formation rates (sSFR = SFR/M) at redshift z m 4 (top left); z ~ 5 (top right); z RJ 6 (lower left) and 
z ~ 7 (lower right). In each instance we use the closest snapshot available to use from our L-Galaxies run of z=3.95, 5.03, 5.82 and 6.97 
respectively. The histogram density plot represents the L-Galaxies galaxy population, with white representing the most dense, and blue 
representing the least. The solid line shows the L-Galaxies median values, and the dashed lines show the 0.16 and 0.84 percentiles. The 
observations are taken from Salmon et al. (2015), the points denotes the median while the error bars reflect the scatter in the observed 
values (not the uncertainty on the median). 


L-Galaxies and Eagle match a similar shape to the obser¬ 
vations, particularly finding good agreement with the slope 
and abundance for low mass galaxies, Illustris predicts 
a slope that steepens with increasing redshift faster than 
what is observed, and over predicts the abundance of low 
mass galaxies at all redshifts. 


5 SPECIFIC STAR FORMATION RATE 

The specific star formation rate (sSFR) is a measure of how 
quickly a galaxy is forming its stars. We present the sSFRs 
at z £ {4, 5, 6, 7} of our galaxy population predicted by our 
model in Figure 6 alongside recent observational measure¬ 
ments from Salmon et al. (2015). We represent the sSFR. of 
individual galaxies by a 2D histogram; the solid line shows 


the median value predicted by our model, averaged over bins 
of 100 or more galaxies. 

The observations are consistent with our model, par¬ 
ticularly for galaxies of mass M ~ 10 9 Mq , across all red- 
shifts. However the observations show a decline in the sSFR 
with increasing galactic stellar mass and we do not iden¬ 
tify the same trend. Instead, all galaxies in our model have 
roughly the same level of activity, regardless of galactic stel¬ 
lar mass. This discrepancy is not surprising: given that the 
models match the observed SFR but under-predict the stel¬ 
lar masses of the largest galaxies, then we would expect this 
result. 

The question remains as to whether the observations or 
the model is at fault, or a combination of both. We could 
boost AGN feedback in the most massive galaxies in the 
model, but this would then reduce the bright end of the 
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Figure 7. Plot to show the evolution of the stellar mass (top left); star formation rate (top right); UV luminosity function (lower left); 
and specific star formation rate (lower right) in the redshift range z = 4 — 7. In each instance we use the closest snapshot available to 
use from our L-Galaxies run of z=3.95, 5.03, 5.82 and 6.97 respectively. 


UVLF. Alternatively, as hinted in the previous section, the 
inferred masses of the highest-mass galaxies may have been 
boosted by observational scatter. At lower masses, there will 
be an observational bias towards the brightest galaxies, and 
so the median sSFR may be over-estimated. 


6 EVOLUTION 

To make it easier to see how the properties of our model 
galaxies change with time, we extract the model predictions 
from each of the four redshifts shown in Figures 2, 3, 5 & 6 
and display them in single panels in Figure 7. 

Concentrating first on the SFR (upper-left panel), we 
see the the knee of the distribution remains relatively un¬ 
changed, at about 20/V 2 Mq yr -1 over this period. However, 
the normalisation of the relation grows and the slope de- 
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creases, such that the number comoving number density 
of galaxies with star-formation rates of O.3 /j _ 2 M0 yU 1 is 
approxiately constant, while that of higher star-formation 
rates in excess of 100/i~ 2 Mq yr -1 grows by several orders of 
magnitude. As might be expected, a similar, but less pro¬ 
nounced, trend is seen in the UVLF, although the knee of 
the distribution is harder to discern. 

In constrast to the SFRDF, the galactic SMDF shows 
only a slight reduction in slope from z = 7 to 2 = 4. Conse¬ 
quently, the comoving number density of low-mass galaxies 
increases by about 1 dex over this time. This is reflected in 
the specfic star-formation rate, which reduces by a factor 
of about 0.5 dex in the same period (as the sSFR is ap¬ 
proximately independent of mass, this conclusion holds for 
individual galaxis, not just the population). 4 


7 CONCLUSIONS 

We have presented the latest liigh-redshift observational 
predictions of the star-formation-rate distribution function 
(SFRDF); UV luminosity function (UVLF); galactic stellar 
mass function (GSMF) and specific star-formation rates (sS- 
FRs) of galaxies from the latest version of the L-Galaxies 
semi-analytic model (Henriques et al. 2014). Our conclu¬ 
sions are as follows: 

(i) We find a good fit to both the shape and normaliza¬ 
tion of the observed SFRDF at z = 4 — 7 (Figure 2), 
apart from a slight under-prediction at the low SFR 
end at 2 = 4, possibly caused by a lack of SFR contri¬ 
bution from merger-driven activity in our model. 

(ii) We find a good fit to the faint number counts for the 
observed UVLF (Figure 3). At brighter magnitudes, 
our predictions lie below the observations, increasingly 
so at higher redshifts. 

(iii) At all redshifts and magnitudes, the raw (unattenu¬ 
ated) number counts for the UVLF lie above the ob¬ 
servations, and so we interpret our under-prediction as 
an over-estimate of the amount of dust in the model 
for the brightest galaxies, especially at high-redshift 
(Figure 4). 

(iv) While the shape of our SMF matches that of the ob¬ 
servations, we lie between the observations at 2 = 4 — 5 
and under-predict at 2 = 6 — 7 (Figure 5). We note, 
however, that both sets of observations are inconsis¬ 
tent with one another, and have, at times, large errors 
attached to them. 

(v) The sSFRs of our model galaxies (Figure 6) show the 
observed trend of increasing normalisation with red- 
shift, but do not reproduce the observed mass depen¬ 
dence, indicating instead that galaxies of all masses the 
same level of activity. It is unclear as to whether this is 
caused by observational bias, or by an under-estimate 
of AGN feedback in the model. 


4 The age of the Universe roughly doubles over this period; thus 
the sSFR measured in terms of this age shows much less variation 
and even at z = 4 is sufficient to double the mass of a galaxy in 
less than a quarter of the age at that time. 


In summary, the L-Galaxies model has mixed success in 
reproducing observations at high redshift. It provides a rea¬ 
sonable match to both the SFR DF and the low-mass end of 
the SMF, but fails to show the observed mass-dependence of 
the sSFR. The predicted UVLF is highly-dependent upon 
an ad-hoc scaling with redshift of the dust model. 

In Yates et al. (2013) we added a detailed model of 
the chemical enrichment in L-Galaxies by adding a de¬ 
layed enrichment from stellar winds and supernovae, as 
well as metallicity-dependent yields and the tracking of 
eleven heavy elements (including O, Mg and Fe). This shows 
promising results in reproducing the mass-metallicity rela¬ 
tion at z = 0, although the chemical enrichment at high- 
redshift remains untested and is something we will look at 
in the future. That then will provide a more realistic predic¬ 
tion of the metallicity of galaxies at early times. In future 
work, we will also add a physically-motivated model for dust 
growth and destruction, and consider the effect of extinction 
from the inter-galactic medium. 
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